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ABSTRACT 

We analyze the clustering of red and blue galaxies from four samples spanning a redshift range of 0.4 < 
z < 2.0 to test the various scenarios by which galaxies evolve onto the red sequence. The data are taken from 
the UKIDSS Ultra Deep Survey, DEEP2, and COMBO-17. The use of clustering allows us to determine what 
fraction of the red sequence is made up of central galaxies and satellite galaxies. At all redshifts, including z = 0, 
the data are consistent with ^ 60% of satellite galaxies being red or quenched, implying that ^ 1/3 of the red 
sequence is comprised of satellite galaxies. More than three-fourths of red satellite galaxies were moved to the 
red sequence after they were accreted onto a larger halo. The constant fraction of satellite galaxies that are red 
yields a quenching time for satellite galaxies that depends on redshift in the same way as halo dynamical times; 
fg ^ (1 In three of the four samples, the data favor a model in which red central galaxies are a random 

sample of all central galaxies; there is no preferred halo mass scale at which galaxies make the transition from 
star-forming to red and dead. The large errors on the fourth sample inhibit any conclusions. Theoretical models 
in which star formation is quenched above a critical halo mass are excluded by these data. A scenario in which 
mergers create red central galaxies imparts a weaker correlation between halo mass and central galaxy color, 
but even the merger scenario creates tension with red galaxy clustering at redshifts above 0.5. These results 
suggest that the mechanism by which central galaxies become red evolves from z = 0.5 to z = 0. 
Subject headings: cosmology: observations — galaxiesxlustering 



1. INTRODUCTION 

In the redshift-zero universe we observe a distinct bimodal- 
ity in the distribution of galaxies; the two populations are 
comprised of blue star-forming objects an d red passively 
evolving objects (Strateva et al. 2001; Blant on et al.l l2003at 
iKauffmann et al.ll200 3b: Madawick et al. 2003). Bimodahty 
is clearly seen at z ~ 1 (Bell et al. 2004; Cooper et al. 20061 
IWillmer et al.l l20()6l) . Recent results detect this bimodality 
out to z = 2 and suggest that the r ed sequence ex ists at even 
higher redshifts (Ki-iek et al. 2008; WilUams et 30 2009. here- 
after W09). How this bimodality was created and has evolved 
is an outstanding problem in galaxy evolution and the his- 
tory of star formation in the universe. In this paper we use 
the clustering of several bimodal galaxy samples spanning 
0.4 < z < 2.0 to test the various physical processes that may 
halt star formation within a galaxy and initiate its transition 
onto the red sequence. 

There are three dominant physical mechanisms by which 
galaxies can migrate from the blue cloud to the red sequence: 
(1) Major mergers can rapidly exhaust the available gas within 
two galaxies in a burst of star formation, resulting in a 
passively evolvin g elliptical galaxy (e.g ., Toomre & Toorrug, 
Il972; Negropo nte & \\Tiitdll983tlMihos & Hernquistlll996r 
[Springel.2000) : t his idea has be e n develo ped in a more sta- 
tistical fashion bv lHopkins et al.l (l2008blah . Since the galaxy 
merger rate is expected to increase with mass, this scenario 
implies a larger red fraction at higher halo mass. (2) Galax- 
ies can be accreted onto a group-sized or cluster-sized dark 
matter halo; in this scenario, the new satellite halo (which 
we will refer to as a subhalo) is no longer able to accrete 
gas, and the gas already contained within the subhalo is sub- 
jected to both tidal forces and ram-pressure stripping. Thus 
star forma tion is predicted to be attenuated in sa t ellite galax- 
ies (e. g.. iGunn & Gotlill972t lAbadi et al.lll999t IWang et all 



effect is obse rved in the color distributions of clusters of 
galax ies (e.g., iButcher & Oemleiill978l 119841: [Hansen et al] 
[2009). (3) There is a halo mass scale above which star 
formation is significantly attenuated, either by virial shocks 
and r e lated processes (Birnboim & Dekel 2003; Kere s et alj 
2005L 12 008!: iDekel & Birnboim .2006; ,Cattaneo etali 12006 



Dekel & Birnboim 2008) or through feedback by active galac- 



120071; Ivan den Bosch et all l2008riKimm et al.1 l2009h . This 



tic nuclei, which is strongly correlated with halo mass 
dCrotoneta l . 2006; Bower etal. 2006; De Lucia & Blaizol 
l2007HSomerville et al.ll2008l) . Many of these studies give rise 
to a critical mass scale of ~ 10'^ /i"' M0, above which shock 
heating within a halo is efficient and gas accretion enters the 
"hot mode" phase. Outstanding questions remain: When do 
these processes begin and how long do they persist? Which 
dominate at high redshift? How is the red sequence built over 
time? 

In this paper, we analyze previously published cluster- 
ing measurements of red (quiescent) and blue (star-forming) 
galaxies for three different surveys at three distinct red- 
shifts. W09 presented clustering results from galaxies in 
the UKIDSS Ultra Deep Survey (UDS). At 1 < z < 2, they 
demonstrate a clear bimodal sequence in the galaxy popula- 
tion between star-forming and quiescent galaxies, with the 
quiescent systems being significantly more clustered. At 
0.7 < z < 1.2, ICoil et alJ (l2008l) presented color-dependent 
clustering measurements from the Deep Extragalactic Evolu- 
tion Probe (DEEP2) spectroscopic survey. Once again, the red 
galaxies showed stronger clustering at all scales than the blue 
galaxies for any mag nitude threshold. U sing the COMBO-17 
photometric survey, Phlep s et"an (12006') (hereafter P06) mea- 
sured the clustering of blue and red galaxies at 0.4 < z < 0.8, 
with similar results. 

Our analysis of this clustering utilizes the halo occupa- 
tion d istributio n (HOD; see, e.g., Peacock & Smith 2000; 
ISeli^HoOd: IScoccimarro et alJl200lt ICoorav & Shethll2002l: 
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Fig. 1 . — Green circles show the number densities of all (red+blue) galax- 
ies in each sample as a function of their mean redshift. The horizontal error 
bars show the redshift range of each sample. The red circles show the num- 
ber density of the red subsample in each survey. The inset panel shows the 
fraction of galaxies that are red in each sample, which, with the exception 
of the bright DEEP2 sample, is roughly 25%. The gray shaded band is the 
number density of galaxies brighter than L* at z = from the SDSS r-band 
luminosity function iBlanton et aU t20Q3bi) ). 



Berlind & Weinberg) l2002l for early works, an d IZhe ng et"an 
20071; Ivan d en Bosch et all 120071: iTinker et al . 2009a for ex- 
amples of more recent implementations of the framework). 
The HOD specifies the connection between galaxies and ha- 
los on a statistical basis, a connection that can be constrained 
through analysis of the two point correlation function ^(r). 

Halo occupation analysis of bicolor clustering provides a 
test of many of the scenarios described above, since each one 
makes a distinct prediction for how red central galaxies oc- 
cupy their halos and how they cluster This enables us to 
address many questions related to the buildup of the red se- 
quence. What fraction of the red sequence is created through 
satellite accretion? What is the quenching timescale for this 
process? At what halo mass scale — if any — do central galax- 
ies transition from the blue cloud to the red sequence? 

Unless otherwise stated, all calculations adopt a flat ACDM 
cosmology c onsistent with the la test constraints from CMB 
anisotropies jDunklev et al.ll2008h . Our cosmological param- 
eter set is (n„„ erg, /J,n.,f^i) = (0.25, 0.8, 0.7, 0.95, 0.045). AH 
distances are in comoving units, and all quoted magnitudes 
assume h= 1. Throughout this paper, we will use the terms 
"red" and "blue" interchangeably with "quiescent/passive" 
and "star forming." 

2. DATA AND ERRORS 

2.1. Clustering in the UDS 

The highest redshift clustering data comes from W09, in 
which the angular clustering, w(9), of star-forming and qui- 
escent galaxies were measured. The galaxies in the UDS are 
selected from observed-frame /T-band photometry, complete 
down to K < 22.4. The sample analyzed here is a photomet- 
ric redshift cut at 1 < z < 2, with a median redshift of z = 1 .4. 
W09 determined the rest-from UVJ colors of the galaxies in 
a UDS fields. Within this color-color diagram, two distinct 
populations are seen out to z = 2, with one population having 
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Fig. 2. — Examples of various central occupation functions for full samples 
and red galaxy subsamples. The thick curve shows a typical (A'cen)M for 
central galaxies that matches the UDS sample. The three thin curves show 
various /rcch models that all have the number of red central galaxies, 20% 
of the number of all central galaxies. The solid curve is a model in which 
/keen = 0.2 at all masses, yielding Mjed/Mtiue = 1- The dotted curve shows a 
model with a broad, smooth transition from halos hosting entirely blue central 
galaxies at low masses to predominantly red central galaxies at the group and 
cluster scale, yielding Mfed/^blue = 3.1. This model is similar to a scenario 
in which central galaxies become red through major mergers. The dashed 
curve shows a model that contains a critical mass scale atM= 10'^ ft-' Mq, 
above which 80% of central galaxies are red, yielding Mred/Mbiue =4.6 . This 
model is similar to a scenario as outhned bv lDekel & Bimboinj (2WB ). 



Star formation rates on average an order of magnitude higher 
than the other W09 found that the clustering amplitude of 
red galaxies is nearly a factor of two higher than that of blue 
galaxies at large scales." 

The errors presented in W09 are bootstrap errors only, ac- 
counting for shot noise but without any estimation of sam- 
ple variance. To place more robust errors on these data we 
use a high-resolution N-body s imulation pe rformed by Mar- 
tin White using his TPM code (White 2002). The simulation 
volume is 720 Mpc per side with 1500^ particles, yielding 
a mass resolution of 7.67 x 10*^ Mq per particle. The cos- 
mology of the simulation is the same as that assumed through- 
out this paper. Using the best-fit HOD model (utilizing W09 
errors only; see in §3 for details on HOD fitting), the halos 
in the simulation were populated with red and blue galaxies 
matching the best-fit occupation functions of each. Although 
the photometric redshift cut in the W08 sample is 1 < z < 2, 
the redshift distribution, A^(z) falls off rapidly after z = 1 .6, and 
the fraction of galaxies at z > 1 .7 is small. The box length it- 
self is not long enough to cover the comoving distance from 
1 < z < 1.7 (roughly 1 /z"'Gpc). To circumvent this prob- 
lem, we reflected the box around the x and y axes, creating a 
volume 4 times larger, then rotated the resulting shape about 
the z axis by an angle of 45°. Within this volume there is 
a cuboid with xyz dimensions 720/V2 x 720\/2 x 720 that 
fully utilizes the original volume of the simulation with no 
double-counting of any structure within the box (here the y 
axis is the line of sight). The extra box length comes at the 
expense of the one other dimension, but the resulting volume 
has sufficient comoving depth to model the W09 data within 
l<z< 1.8. 
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Fig. 3. — Panel (a): Clustering measurements from lWilliams et al] j2009l) for quiescent and star-forming galaxies in the Ultra Deep Survey of UKIDSS within 
the redshift range I < z <2. Circles represent the quiescent subsample, while squares represent the star-forming galaxies. The solid and dotted curves show the 
best-fit HOD model for quiescent and star-forming galaxies, respectively. Panel (b): Halo occupation functions from the best-fit model. The thin solid curves 
indicate the HOD for all galaxies in the sample; the horizontal curve with a Gaussian cutoff at A/ < lO"^ ft"' Mq represents all central galaxies, while the 
power-law curve represents satellite galaxies. The thick, red curves indicate the fraction of these galaxies that are quiescent (red). The dashed curves represent a 
model in which /Rcen = 0.22 independent of halo mass. This model is within the la lower limit on the HOD from the clustering data. 



Given the angular area of the the UDS, 0.7 deg^, 140 in- 
dependent lines of sight could be obtained from the simu- 
lation. In each hne o f sight, we calculate w(9) with the 
iLandv & Szalavl(ll993 l) estimator, as in done in W09. Ideally, 
a full error analysis should include the entire covariance ma- 
trix of both the red and blue clustering measurements, which, 
since they are measured from the same patch of sky, will be 
correlated. The 140 realizations are not sufficient to obtain 
a reliable estimate of the covariance, so throughout this paper 
we use diagonal errors only on the W08 data. This renders the 
parameter constraints obtained somewhat suspect. To aide in 
the constraining power of these data, and to help make up for 
the lack of an estimate of the covariance between the red and 
blue galaxies, we also use the relative bias of red and blue 
clustering. This quantity is less effected by sample variance 
because the amplitudes of red and blue galaxies in a given 
patch of sky will move in concert if the o verall clustering o f 
that patch varies from the universal mean (Seljak et al. 2009). 
Taking the ratio of w{d) at each of the four data points at 
9 > 100", where the data are firmly in the two-halo regime, 
the relative bias bi-ed/buue = 1.42. Using the same approach 
on w(9) from the mocks, the variance of /jred/^Wue is 11%. We 
also incorporate the error in the number density of each sub- 
sample, which we estimate from the variance of the mocks. 
These errors are 12% for the red galaxies and 14% for the 
blue galaxies. 

Finally, to model w(9), information on the redshift distribu- 
tion, N{z), is required. The photometric redshift distribution is 
close to a top-hat function from 1 < z < 1 .8, but due to uncer- 
tainties in the photometric redshifts, the true N(z) is somewhat 
broader An estimate of the underlying N{z) using the errors 
in the photo-z's, kindly provided by R. Williams, is used in 
all analytic calculations. The methodology for obtaining the 



estimated N(z) is described in W09 and lOuadri et al.l ( 120081) . 

2.2. Clustering Measurements from the DEEP2 Survey 

Measurements of the clustering of blue and red galax- 
ies a^j_£;;^. 8 inthe DEEP2 survey have been published 
bv Co il et al.l (l2008h . Being a spectroscopic redshift survey, 
DEEP2 is superior to the UDS in that the true redshift-space 
correlation function can be estimated. This quantity, referred 
to as C(r^,7r), where is the projected separation of galaxy 
pairs and tt is the line-of-sight separation, can be integrated 
along vr to obtain a measurement of the projected galaxy cor- 
relation function, Wp{rp). The use of a projected quantity 
ameliorates the effect of redshift-space distortions on both 
large scales from coherent infall into overdense regions (e.g.. 
Kaiser 1987; Fisher 1995; Hamilton 1998; Scoccimarro 2004]) 
and small scales from non-lin earities a nd vir ial motions of 
satellite galaxies (e.g.. lDavis & P eebles . 1983h . The quantity 
Wp(rp) is only free of redshift space effects if the projection 
is over a sufficiently large line-of-sight distance (formally, 
the projection of the real-space and redshift-space correlation 
functions are identical if the projection is extended to infinity). 
Due to the small sample size of DEEP2, 3 square degrees, the 
statistics are not goo d enough to estim ate £,irp,TT) out to large 
values of tt, thus the Coil et al.l (12008 ) measurements are inte- 
grated over the range ±20 Mpc. This tt^^^ is large enough 
to eliminate non-linear redshift-space effects, but large-scale 
infall will contribute to the measured amplitude of Wp(rp) at 
rp > 5 /i"'Mpc. 

Our goal in this paper is to probe the galaxy and halo mass 
scale at which the transition to the red sequence begins. Thus 
we focus on the DEEP2 results for the magnitude threshold 
sample of Mb < -19.5, the faintest sample published. For the 
blue galaxies, the median redshift is 0.81, while for the red 
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galaxies the median redshift is 0.77. In practice we use z = 
0.79 as the redshift of the combined red-blue analytic model, 
noting that our results depend little on the assumed redshift 
in the range specified. As a consistency check on our results, 
we also consider a bright sample of galaxies with Mb < -20.5. 
The median redshifts for the blue and red subsamples are 0.99 
a nd 0.88, res p ective ly. 

ICoil et ai] (l2008l) estimated errors from the sample-to- 
sample variance from 10 independent fields. Since the vol- 
ume of this sample is ~ (100/i~' Mpc)-', sample variance is a 
concern. As with the UDS sample, we also include the uncer- 
tainties in the number densities for each sample. We estimate 
these uncertainties from mocks created from the sample N- 
body simulation as above, but now each mock is a cube with 
the same volume as the observational sample. The mocks are 
created from using the best-fit HOD parameters under the as- 
sumption of no errors in the number density. The uncertain- 
ties in the abundances are 11% and 13% for the blue and red 
subsamples, respecti vely. 

All the samples in Coil et al.l (1200 8^) are separated into blue 
and red subsamples by a tilted cut in color-magnitude space, 
(U-B) = -0.023(Mb + 2 1.62)+ 1.035. Using this cut, - 27% 
of the galaxies in the faint sample are red, while ^ 36% of 
the objects in the bright samples are classified as red. As op- 
posed to the UVJ method of defining samples in the UDS, a 
single color cut can be subject to dust contamination. We will 
discuss the possible systematics of this effect in §6. 



2.3. Clustering Measurements from COMBO-17 

The COMBO-17 survey is a photometric survey covering 
0.78 deg^ with 17 broad- and medium-ban d filters, produc- 
ing high-precision photometric redshifts (see Wolf et al.l2003L 
|2004 for details of the survey). Galaxies are identified in 
observed-frame /?-band, from which P06 created a volume- 
limited sample of galaxies down to Mb < -18. The photo- 
metric redshifts are accurate enough to yield measurements 
of £,irp,TT) and Wp(rp) at a redshift of z = 0.6, presented in 
P06. Unlike the DEEP2 sample, in which Wp{rp) is obtained 
from integrating ^(r^, tt) out to 20 Mpc, P06 integrate out 
to 100 Mpc using the direct measurements, then continue 
out to 200 Mpc using a linear theory model set to match 
the amplitude of the data at tt = 100 /!"'Mpc. With this ap- 
proach, redshift-space distortions are eliminated in the esti- 
mate of Wp(rp), provided that the linear theory model is accu- 
rate. Error bars are estimated from jackknife sampling of the 
COMBO-17 survey. Although P06 present measurements of 
the covariance matrix of both the red and blue subsamples, the 
matrices themselves are too noisy to be properly inverted. We 
use the diagonal errors obtained from the jackknife method in 
our analysis, but it should be noted that the matrices presented 
in Phleps et al. (2006) show that Wp(rp) becomes correlated 
outside of rp ^ 1 Mpc. We assume the same uncertainties 
in the number densities as for the faint DEEP2 sample. 

Utilizing the wide wave l ength coverage of the 17 photo- 
metric bands, Phlep s et"aLl (l2006l) used spectral-energy dis- 
tribution (SED) fitting to break their sample of galaxies into 
active and passive types. The template spectra used in the fit- 
ting comprise a two-dimensional age/reddening sequence to 
remove the effect of dust on the classificatio n of galaxies. A 
U — V color cut based on the prescription of Bell et al] (l2004l) 
is then used to divide the sample into red and blue galaxies. 
This cut is a function of My and redshift. 



2.4. Comparison of Surveys 

These three surveys make up a heterogeneous collection of 
samples that vary in the bands, magnitude limits, and def- 
initions of red and blue. It is perhaps more productive to 
look at the number densities rather than the luminosities — the 
number density of a threshold sample of galaxies correlates 
strongly with the mean halo mass scale being probed. The 
halo mass function also evolves from z = 1.4 to z = 0.6; the 
abundance of massive halos increases dramatically, while the 
number of M < lO" halos evolves only weakly. At the red- 
shift of the UDS, the nonlinear halo mass scales — where the 
bias is approximately unity — is 1.7 x 10"^ H'^Mq, while at 
z = 0.6 it is 3.1 X 10" /;"' Mq. At all redshifts, the halo mass 
scale probed by each sample is above this nonlinear scale, 
implying that any difference in the halo masses probed by red 
and blue galaxies will produce differences in the large-scale 
clustering. 

The mean redshifts, redshift ranges, and number densities 
of the galaxy samples are shown in Figure [T] Ignoring the 
bright DEEP2 sample for the moment, the number density in 
our samples increases by roughly a factor of 3 from z 1 .4 to 
z ^ 0.6. Thus COMBO-17 is probing the clustering of lower- 
mass halos than either DEEP2 or the UDS. However, the frac- 
tion of galaxies that are red or passive is roughly constant at 
- 25%. 

3. METHODS 

3.1. The Halo Occupation Distribution for a Threshold 
Sample of Galaxies 

Before describing our model for the occupation of galaxy 
samples defined by color, we start by defining the halo occu- 
pation of all galaxies, red and blue together For galaxy sam- 
ples that are complete down to a given luminosity threshold, 
halo occupation is broken into two distinct parts: galaxies that 
reside in the center of a dark matter halo and satellite galaxies 
that are within the virial radius of the halo but are distributed 
throughout the halo. 

For central galaxies, we parameterize the central occupation 
function as 



1 

{Ncen)M - 2 



^^^^ /logM-logM, 



(1) 



Equation ([T]i yields a smooth transition between halos that are 
too small to contain galaxies bright enough to be included 
in the given sample (M ^ M^i„) and halos that are massive 
enough such that they will always contain a galaxy at their 
center bright enough to be included (M ^ M^i„). Formally, 
Mmin is defined as the halo mass at which a galaxy has a 50% 
probability of containing a central galaxy in the sample. The 
parameter uiogM controls how rapid the transition is between 
zero and one central galaxies. Physically, aiogM represents the 
scatter in halo mass at the luminosity threshold defined by the 
sample, under the assumption that this scatter takes the form 
of a lognormal distribution. Equation ([T]i is the mean number 
of central galaxies; because there can only be one or zero cen- 
tral galaxies in a halo, the scatter about that mean is defined by 
a nearest-integer distribution. Low-redshift estimates of the 
mass-luminosity scatter yield values between 0.2 and .6, de- 
pending on lum inosity (Tinker et al. 2006, 2007; Zheng^etaD 
I2007t Ivan d en Bosch et al. 2007t iMore et al.ll2009 ). For the 
UDS and COMBO-17 samples we set uiogM = 0.3 while for 
both DEEP2 samples we set tiiogM = 0.6. The motivation 
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Fig. 4. — Panel(a): Measurements of the projected correlation function for galaxies with Mb < -20.5 in DEEP2 bv lCoil et alj<200^ . The sample has a median 
redshift of z ~ 0.9. The circles and squares represent the measurements for red and blue galaxies, respectively. The solid and dashed curve represent the HOD 
model fits for red and blue galaxies, respectively. Panel (b): The best-fit occupation functions for all galaxies in the sample (thin solid curves) and red galaxies 
(thick curves). The model results imply that there is no threshold mass scale at which central galaxies preferentially become red; the higher clustering amplitude 
of red galaxies is produced by the high fraction of red satellite galaxies. 



for the higher a\ogM for DEEP2 springs from the observed 
clustering amplitude of DEEP2 galaxies; the large-scale am- 
plitude of Wp(rp) is somewhat lower than predictions using 
the WMAP5 cosmol ogy and a reasonable bias model (see 
IWetzel & Whit3l2009,) . This could imply lai- ger magnitude 
errors or sample variance in the data. Allowing (TiogM to be a 
free parameter drives the best-fit values to physically unreal- 
istic regions on parameter space. Fixing these values, rather 
than allowing them to be free parameters, does not bias our 
results. 

For satellite galaxies, we parameterize the mean satellite 
occupation function as 



(A^sat)M = 2 



1+erf 



logM-log2M„ 

TlogM 



M 



(2) 



At halo masses well above the minimum mass scale, the 
number of satellite galaxies scales as a power-law with in- 
dex Qfsat and normalization Msat- This is well-motivated 
from the results of hydrodynamical cosmological simula- 
tions (e.g.. White et al. 2001; Berlind et al. 2003; Zheng et al, 
20051) as well a high-r e solution coUis i onless simulations (e.g.. 



Kravtsov et all l2004t IConrov et al.1 l2006t IWetzel & Whita 
20091) . In equation ©, the satellite occupation function has 



a cutoff of the same functional form as (A'cen)M, but with a 
transition mass a factor of two larger that for central galaxies. 
This cutoff prevents halos with a low probability of containing 
a central galaxy from having a higher probability of contain- 
ing a satellite galaxy in the sample, and is also motivated from 
the numerical results listed above. We assu me that the scat- 
ter about the mea n (A^sat)M is Poisson (Krav tsov et al.ll2004t 
IZheng et al.ll2005l) . 

We set asat = 1 for all samples. There is strong theoreti- 
cal and observational evidence for such a prior High resolu- 



tion dark matter simulations find that the subhalo mass func- 
tion is nearly self-similar dGao et al.l 120041; iDe Lucia et alj 
2004; Kravtsov et al. 2004|), thus if satellite luminosity is 
tightly correlated with subhalo mass' the number of satel- 
lites above a given luminosity should scale linearly. Both cos- 
mological hydrodynamic simulations and semi-analytic mod- 
els of galaxy formation find that this scaling is not broken 
with the inclusion of ba ryonic physics (White et al. 200]]; 
IZheng et all l2005t fCroton et all I2006L iHarker et al.1 12007). 
Observation ally, there is som e discrepancy in the measured 
in clusters. iLin et al.l (120041) fin a value of asat < 1, while 
iKochanek etaP (l2003h find a^at = 1. 10 ± 0.09 (for a sample 
that contains significant overlap with the Lin et al. 2004 sam- 
ple). IColhster & Lahavl (l2005l) . using a much larger sample of 
galaxy gro ups from the 2PIG G catalog (Eke et al..i2004 ). find 
ttsat = 0.99. lYang et akl (|2008|) find asat ~ 1 from a group cata- 
log constructed from the spectroscopic galaxy sample in DR4 
of the SDSS. In the largest sample of galaxy clusters to date, 
the maxBCG cluster catalog (Koester et al. 2007), the results 
are also consistent with asat = 1 (Hansen et al. 2009) . We will 
discuss the effect of leaving asat as a free parameter for both 
red and blue samples in §6. 

3.2. Separating Galaxies into Red and Blue 

Many previous studies have presented halo occupation 
modeling of red and blue galaxy clustering ( Scranton 2003]; 
'Magliocchetti & Porciani 2003; Zehavi et alj l2005i: P06; 
Tinker et al. 2008b; Skibba & Sheth 2009). We will present a 
somewhat modified approa ch to color-dependen t HOD mod 
eling that is an extension 



aach to color-de penden t 
of [Tinker et all ( (2009al) . 



' As not ed in 'Nagai & Kravtsov' (TOOf), 'Vale & Ostriker' 0003) . 
IConrov et all (2006), Wans et al. (2006). and Wetzel & White (2009), satel- 
lite luminosity or mass correlates best with the mass (or circular velocity) of 
the subhalo at the time it is accreted, rather than the time it is observed. 



6 




[h-' Mpc] log M, [h-' MJ 

Fig. 5. — Panel (a): Measurements of the projected correlation function for galaxies with Mb < -19.5 in DEEP2 by ICoil et al J 120081) . The sample has a 
median redshift of z ~ 0.8. The circles and squares represent the measurements for red and blue galaxies, respectively. The solid and dashed curve represent 
the HOD model fits for red and blue galaxies, respectively. Panel (b): The best-fit occupation functions for all galaxies in the sample (thin solid curves) and red 
galaxies (thick curves). The asymptotic fraction of central galaxies that are red, /Rmax, is set to be 0.25 to match the results from Figure|4] As with the bright 
sample from DEEP2, the HOD results for this sample imply that there is no preferred mass scale for red central galaxies; increased clustering of red galaxies is 
due mostly to /Rsat , which is nearly twice as high as the overall fraction of red galaxies in this sample. 



For satellite galaxies, we utilize a simple model to sepa- 
rate red and blue subsamples: a constant fraction of satel- 
lites, /Rsat, are red, independent of halo mass . Previous 
models for low -redshift color-dependent data, ie Zehavi et al.l 
(|2005), allow this red satellite fractio n to vary with halo m ass. 
But in most samples investigated bv lZehavi et al.l (l2005h and 
[Tinker et alj ([2008b). as well as our own tests, the best-fit 
models have little to no variati on of /Rsat with M . In the 
DEEP2 galaxy group catalog of iGerke et al.l ( !2007j) . the red 
fraction is independent of group richness, supporting our as- 
sumption that the /{sat is independent of halo mass. We will 
discuss the possible biases of this assumption in §6. 

For central galaxies, we implement a model with con- 
siderably more flexibility. Because the the efficiency of 
star formation quenching for central galaxies may be halo 
mass dependent, the minimum mass scale for red central 
galaxies may be much higher than that for blue galaxies 
in a luminosity-threshold sample. Current models predict 
that this transition could be sharp (e.g., Croton et al. 200^ 
iDekel & Birnboiml 120061: Cattane o et al.1 120061). or fliat it 
could be quite broad (e.g.. .Bower et al.ll2006l: iHopkins et al.l 
I2008al: ISomerviUe et al.ll2008h . At high redshift. massive ha- 
los may be more efficient at forming stars than in the present 
universe, thus a hig h fraction of high-mass halos may be blue 
jDekel et al.l l2008l) . We parameterize all these possibilities 
with a red central fraction of the form 



/RcenW = /Rmaxexp 



-(3kM„ 



M-(3M„ 



(3) 



when M > jSM^^^, and /Rcen = at lower masses. The pa- 
rameter K governs how sharp the transition is between halos 
hosting no central red galaxies and halos having an asymp- 
totic probability of /R^ax of having a red central galaxy. The 



parameter [3 sets the overall shift in the mass scale for red cen- 
tral galaxies by setting the mass scale below which no halos 
host red central galaxies. 

Figure |2] shows several examples of the red occupation 
function, defined as /Rcen x (^cen)M- The thick solid line 
shows an example of {Ncen)M, the central occupation function 
for all galaxies, for the UDS sample. Mmin for this sample is 
roughly lO"'^ /i~'Mq. The thin solid curve shows a model 
in which (/Rmax,K,/3) = (0.2,0,0.1), which yields a constant 
/Rcen = 0.2 at all halo masses. A quantity of interest is the 
ratio of the average halo mass of a red central galaxy to the 
average halo mass of a blue central galaxy, defined as 



^red ^ «Rcen/ dMMn{M)f^^^,, X (A^cen)M 
Mblue "Been / dMMn(M)( 1 - /Rcen) X (A^cen)*? 



(4) 



where n(M) is the halo mass function, and WRcen and necen are 
the number density of red and blue central galaxies, respec- 
tively, that are calculated from the HOD by, eg. 



"Rcen = y t/Mn(M)/Rcen X (A^c, 



(5) 



For the thin solid curve, Mied/Mbiue = 1 by definition. Red 
central galaxies are essentially a random subset of all cen- 
tral galaxies. The dotted curve shows a model in which 
(/Rmax,K,/3) = (0.75,24,0.75). This model produces a very 
broad transition between halos having a negligible probabil- 
ity of containing a red central galaxy to central galaxies being 
majority red. The number density of red centrals is the same 
as the solid curve, but now the mass ratio between red and 
blue central galaxies is 3.1. This function would represent a 
merger-induced red central model. The dashed curve shows a 
model in which (/Rmax, k,/3) = (0.8,0.057,2.5). In this model. 
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the transition between blue and red central galaxies is nearly 
instantaneous at M = 2.5Mmin- Once again, WRcen is the same 
as in the previous two models, but M,ed/Mbiue is increased to 
4.6. A function of this form represents the critical mass scale 
scenario. 

3.3. Calculating Observables with the Model 

To calculate the galaxy autocorrelation function, ^(r), from 
a given H OD, we use th e analytic model described in the Ap- 
pendix of [Tinker et alj ([2005) (see also Zheng 2004). As a 
brief description, the correlation function is broken into two 
parts: a one-halo term in which pairs of galaxies reside within 
a single halo, and a two-halo term in which pairs come from 
two distinct halos. The one-halo term dominates the corre- 
lation function at separations r < 1 /z"' Mpc, while the two- 
halo term contributes nearly all pairs at r > 1 /;"' Mpc. The 
shape of the one-halo term is influenced by the radial distribu- 
tion of satellite galaxies, which we assume to be the same as 
the dark matter S pecifically, we use the halo concentration- 
mass relation from lZhao et al.l (l2008l) . More importantly, the 
one-halo term is sensitive to the overa ll fraction of galax ies 
that are satelhtes (se e the discussions in lZheng et al.l200^ and 
[Tinker et alj|200 9a). The satellite fraction is influenced both 
by the shape of the mass function and the form of the HOD. 
More satellites will increase the large-scale bias of a sample, 
but because the number of satellite-satellite pairs within a sin- 
gle halo increase as (A^sat)M, the effect of the small-scale clus- 
tering is more dramatic. At large scales, the shape of the two 
halo term is the same as the linear matter correlation func- 
tion. At intermediate scales, r < 3 Mpc, halo bias deviates 
from a simple li near approximation, and a scale-dependent 
term is required dTinker et al.ll2009"bl) . When implementing 
the analytic model for ^(r), we use the halo mass function of 
[Tin ker et al. (2008a) and the halo bias function of Tinker et al. 
( l2009b) . We assume that a halo is defined to be a spherical 
object with a mean interior density 200 times the background 
density. 

For the UDS, the angular correlation function, w{0), is mea- 
sured. This quantity is defined as 



w{e) = j dzN^(z) -^j dx^ [Vx^ + r^e^j , (6) 

where N{z) is the normalized redshift distribution of the 
galaxy sample, r is the comoving radial distance at redshift 
z and dr/dz = (c / Hp) / \/Q,„(l + z Y + As sta t ed in §2.1, 
we use the estimate of N{z) from IWilliams et ail ( l2009h that 
incorporates photometric redshift errors. 
For DEEP2, we calculate the projected correlation function, 

Wp(rp) = 2 / ^(rp,'K)d7r. (7) 
Jo 

As mentioned in §2.2 we set TTmax = 20 /i"' Mpc, the same as 
in the measurements. For the one-halo term, 20 /;"' Mpc is 
sufficient to eliminate redshift-space effects, therefore we use 
the isotropic real-space one-halo term in equation (Q- For the 
two-halo term, coherent infall is still a concern. We t ake this 
effect into account by using the linear theory model of lKaisej 
(Il987h when calculating the two-halo term. Although linear 
theory is not a fully descriptive model of ^(rp,7r) at separa- 
tions of a few /i"' Mpc, for a projected quantity like Wp{rp) 
we find that it compares well to galaxy mocks constructed 
from N-body simulations in which Wp(rp) is estimated with 



TTmax = 20 enforced. See i Hamiltonl (Il998h and iHawkins et ak 
(2003) for a thorough discussion and application of the lKaisei 
( 1987) model to an unparameterized correlation function. 

Because the COMBO- 17 data integrate to TTmax = 200 
/i"' Mpc, the effect of redshift-space distortion is negligible. 
To calculate Wp{rp) for the COMBO-17 data, no accounting 
of redshift-space effects are required and we use the isotropic, 
real-space ^(r) in equation |7] and integrate along the line of 
site until convergence. 

For each sample, we use the Monte Carlo Markov Chain 
(MCMC) method to probe the likelihood distribution of the 
full parameter space. We have six free parameters. The first 
two, Mmin and Msat, specify (A^) m for the full sample of galax- 
ies. We have one parameter, /Rsat, to split the satellite galaxies 
into red and blue subsamples. The last three parameters, /3, k, 
and /Rmax, determine the fraction of central galaxies that are 
red as a function of halo mass. These parameters specify not 
just ^(r) but also the abundance of blue and red galaxies. We 
minimize the total of the model, which we obtain from the 
sum of the for the red and blue clustering measurements 
and the for the galaxy abundances. The best-fit parameters 
for each sample are given in Table 1 . 

4. RESULTS 
4. 1 . Clustering in the UDS 

Figure [3^ shows the best-fit model against the 1^(6*) data, 
and Figure [3}' shows the HOD for all galaxies and the quies- 
cent subsample. The relative bias between the star-forming 
and quiescent galaxies is high enough that the best-fitting 
model is on the extreme edge of allowed parameter space; 
to increase the clustering of quiescent galaxies relative to the 
star-forming subsample, the best-fit /rcbh rapidly rises from 
to 1 at a "critical" mass threshold of M = 10'^^ /i"' Mq, yield- 
ing a halo mass ratio of red to blue centrals of Mred/A^biue ^ 6. 
Because of the high central mass ratio, the clustering ampli- 
tude of red galaxies is boosted significantly. This fact, along 
with the high fraction of centrals that are red, forces the frac- 
tion of red satellite galaxies down to /Rsat ~ 22%, which is 
lower than the overall red fraction of galaxies. 

However, there is a strong degeneracy axis between 
M,ed/Mbiue and /Rsat- Bccausc the number of satellites in- 
creases linearly with halo mass, increasing /Rsat also increases 
the large-scale bias of red galaxies without any increase the 
halo mass scale of red central galaxies. This, combined with 
the large error bars on the data, make the constraints on the pa- 
rameters of the HOD model weak. Within the 2a confidence 
limit, a model in which /Rsat = 0.59 and Mi-ed/Mbiue = 1-01 is 
acceptable. This model is shown with the dashed line in Fig- 
ure |3]3. We will discuss this degeneracy and its implications 
in more detail §5. 

4.2. Bright galaxies in DEEP2 

Figure |4] presents results for DEEP2 galaxies in the Mb < 
-20.5 sample. Figure |4^ shows both the measured cluster- 
ing signal and the best-fit HOD model. Figure ^ shows the 
HOD for both the full galaxy sample and the red subsample. 
At all scales, the clustering of red galaxies is higher relative 
to blue galaxies. However, the Wp(red)/wp(blue) ratio is sig- 
nificantly small er than tha t found in lower redshift samples 
such as SDSS (IZehavi et al. 2005; Li et al. 2006) and 2dF- 
GRS (Norb erg et al.l2 002; Madgwick et al. 2003). In fact, the 
the increased clustering of red galaxies in the two-halo term 
is best accounted for by a model in which red central galaxies 
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Fig. 6. — Panel (a): Measurements of the projected correlation function for galaxies in COMBO-17 bv lPhleps et al.lf2006l) . The median redshift of the sample 
is z ~ 0.6. The circles and squai'es represent the measurements for red and blue galaxies, respectively. The solid and dashed curve represent the HOD model 
fits for red and blue galaxies, respectively. Panel (b): The best-fit occupation functions for all galaxies in the sample (thin solid curves) and red galaxies (thick 
curves). The halo mass ratio of red central galaxies, Miej/Mjjiue, is 1.27, in agreement with the results from both DEEP2 samples. 



have no dependence on halo mass (Mred/Mbiue = 104), but 
with a high fraction of satellite galaxies being red (/Rsat = 
0.66), twice the global fraction of red galaxies in this sam- 
ple. Unlike the analysis of the UDS clustering, the constraints 
on the model parameters are strong. 

The value, at 64.4 for 29 degrees of freedom, is not ideal 
but we note that without a proper covariance matrix the ^ 
values are not fully robust. The for the blue galaxies (39. 1 ) 
is significantly larger than for the red (18.3). Visually, it ap- 
pears that this is a result of the poor fit to the blue WpiXp) in the 
one-halo term, but point-to-point scatter is also a concern; two 
points, the second and fourteenth, contribute a x^ of 22. To 
improve the fit at small scales, the amplitude of the one-halo 
term for both red galaxies and blue galaxies would need to in- 
crease. This is obtained by increasing the number of satellite 
galaxies, but it would come at the cost of increasing the am- 
plitude of the large-scale bias as well. The best-fit model is 
somewhat above the data at large scales (vp > 2 /i"' Mpc), and 
increasi ng this disparity would decrease the quality of the fit 
(see also lWetzel & Whitell2009h . IZheng et al.l (l2007l) . also us- 
ing HOD modeling, produced fits that accurately reproduced 
the large-scale clustering of DEEP2 galaxies, but this result is 
due to the choice of cosmology: the WMAPl-type cosmol- 
ogy utilized in Zheng et al. (2007) has significantly less large 
scale power due its lack of tilt and hi gher matter d e nsity. We 
find that we are able to reproduce the IZheng et all (l2007h re- 
sult when implementing the WMAPl cosmology. 

In this paper we are most interested in probing the halo 
mass scale at which star formation for red central galaxies 
becomes quenched, if there is one. Analysis of bright galax- 
ies probes a higher halo mass scale, but this sample is useful 
for a number of reasons. First, it is a consistency check on the 
results obtained from the faint galaxies. Second, because the 
samples in Coil et al, (2008J are threshold samples, we know 
that the /Rcen for the faint sample can never be lower than /Rcen 



for the bright sample at a given halo mass. It is also highly 
unlikely that /rcbh for the fainter sample would be higher than 
that of the brighter sample at M > lO'** Mq unless magni- 
tude errors were extreme. The bright sample is fully contained 
within the faint sample; the central galaxies in high-mass ha- 
los are the same galaxies that produce the clustering in the 
bright sample. Thus we fix the value of /Rmax in the faint 
sample to be the the best-fit value of from the Mb < -20.5 
sample, 0.25, ensuring that /rceii for the faint sample will be 
the same as that in the bright sample in the overlapping halo 
mass range. 

4.3. Faint galaxies in DEEP2 

Figure|5h presents the clustering of Mg < -19.5 in DEEP2, 
along with the model with the lowest x^ from the MCMC 
chain. As with the bright galaxies, the clustering of red galax- 
ies in this sample is stronger than that of the blue galaxies at 
all scales (the last three data points in the red Wp{rp) that have 
low amplitudes are most likely a result of sample variance), 
but the relative bias between red and blue samples a t large 
scales is only moderately above unity dCoil et alj|2008l calcu- 
late a relative bias of 1 .25). Thus, the best-fitting HOD in Fig- 
ure [Sj^ is similar to that of the bright galaxies; /rcbh is nearly 
independent of halo mass, yielding Mjed/Mbiue = 1-48. The 
central occupation function for red galaxies contains a sharp 
cutoff at M= lO" "* /!~'M0, which is the feature that causes 
the value of Mred/A/biue to be larger than one. This cutoff is 
required by the constraints on the number density rather than 
on the clustering; without such a cutoff, the space density of 
red galaxies would be too high. The fraction of satellites that 
are red is 52%, twice the overall red fraction of this sample. 

Visually, the fit for this sample appears to match the data 
better than with the brighter sample, but the x^ is roughly 
the same, at 63.7 (although for one less degree of freedom 
because /Rmax has now been fixed to 0.25). The larger num- 
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Fig. 7. — Constraints in the A/red/'Wbluc-/Rsat plane for tfie four samples analyzed in this paper. Contours represent 1- and 2 — tr constraints. In panel (a), 
showing the constraints from the UDS clustering data, the degeneracy between the central mass ratio and the fraction of red satellites is clear. For the other three 
samples, Mred/Mbiuc is consistent with unity, with a red satellites fraction of ~ 0.6. In Panel (c), showing the constraints for the faint DEEP2 sample, the circle 
with errorbars represents the results of a model in which /Rmax is allowed to be a free parameter. In the contours, /Rmax = 0.25 to agree with the results from the 
bright DEEP2 sample. EiTorbars are 2cr. 



ber density of this sample yields smaller error bars in the 1- 
halo term, and once again point-to-point variance is larger 
than the errors for several data points: two data points ac- 
count for 27% of the total x^- When using the WMAPl cos- 
mology to model the red and blue, the of the fit is 28.9 
(yielding x^/ v = 0.93), a significant improvement on the re- 
sults discussed above. However, the constraints on /Rsat and 
Mred/Mbiue are unchanged. An additional difference between 
the IZheng et al. (2007) analysis and that of this paper is the 
freedom in Osat- We will discuss this is more detail in §6. 

As a check on these results, we have performed the HOD 
analysis of the faint sample without the prior on /Rmax ob- 
tained from the analysis of the bright sample. The results are 
consistent within their 1-cr contours, which we will show in 
Figure |7] In this test, /rs^i was nearly unchanged, while the 
best-fit Mi-ed/Mbiue decreased to be 1.0, the same as for the 



bright galaxies. The best-fit value of /Rmax is 0.16, rather than 
0.25, which removes the need for the sharp cutoff in /rcbh at 
10" "* Mq in order to match the number density. 

4.4. Clustering Analysis for COMBO- 17 

Figure|6t and|6j3 show the clustering results from COMBO- 
17 and the resulting halo occupation functions, respectively. 
The Wp{rp) measurements for red and blue galaxies are similar 
to those in DEEP2; at small scales, the clustering of red galax- 
ies is strongly enhanced relative to blue galaxies, but in the 
two-halo term the difference in the clustering is minimal, with 
sample variance muddling the comparison at rp > 10 Mpc. 
This results in a best-fit HOD that is the same as those found in 
the DEEP2 sample; M^ed/Muue = 1.01 and /Rsat = 0.58, which 
is more than twice as high as the overall red fraction of the 
sample. 

The reduced x^ value of the best fit model is 1.9. The x^ for 
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Fig. 8 . — The subhalo mass function from a numerical simulation pre- 
sented in IWetz el & Wh it3 (2009). Here Ms„(, is the mass of the subhalo at 
the time of accretion. The thin curves represent the N-bo dy results, while the 
thick curve is the fitting function given in Equation (TO). For the numerical 
results, each color represents a different host halo mass bin, ranging from 
logMhost = 12.4 to logMhost = 13.9. For each bin in Mhost. the six lines show 
results for six redshift outputs ranging from 7, = 0.8 to z = 1.6. The turnover 
at low Msiib/Mhost for each mass bin is due to the fixed particle limit in the 
subhalo catalog; we use only subhalos with more than 2000 particles at in- 
fall to avoid artificial numerical disruption. The overlapping results for each 
mass bin show that the subhalo mass function is universal over a wide range 
of halo mass and redshift. 
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Fig. 9. — The distribution of subhalo ages, defined as the time since a 
subhalo was accreted onto its parent halo. Results are from the simulations 
of Wetzel & White ( 2009). The y-axis indicates the fraction of subhalos with 
accretion times longer than the value t. As redshift increases, the mean age 
of a subhalo becomes younger because dynamical friction is more efficient 
at higher redshifts due to the higher density of most halos. The halo catalogs 
span the redshift range z = to z = 1.6, but in this plot we show only four 
output redshifts to avoid crowding. 



the blue Wpijp) is 14.7 for 25 data points. The source of the 
high is localized to the transition region between 1-halo 
and 2-halo clustering in the red Wp(rp); half of the total 
value is accrued between 0.5 < < 1.1, where the clustering 
amplitude drops rapidly. P06 have also performed HOD anal- 



ysis of these data, obtaining a much better ^ value than that 
found here, 35.9 versus 88.3. There are several reasons for 
the lower value from the P06 analysis. First, P06 have a 
different cosmological model with no tilt, which reduces the 
large-scale power. Second, P06 leave erg as a free parame- 
ter, increasing the freedom of their cosmological model. This 
parameter is left free for both the red and blue subsamples, 
yielding best-fit values of a% of 0.84 ± 0.08 and 1.19 ± 0.09 
for the red and blue subsamples, respectively. The lower (Tg 
value for the red model fit is required by the low relative bias 
of blue and red galaxies in the data. Third, and most impor- 
tant, P06 implement a physically unrealistic HOD model. P06 
assume the same form of the HOD for both blue and red sub- 
samples: a power-law in which (A^) m = (M/Mq)'^ for M > Mq. 
The power-law index (3 is left as a free parameter. This model 
is unrealistic in that it does not reproduce known aspects of 
halo occupation parameterized in equations ([T]i and a 
self-consistent separation of central and satellite galaxies; a 
"shoulder" in (A^) m in the halo mass range between Mmin and 
Msat (nearly a factor of 20 in halo mass), in which the scat- 
ter in (A^) M is sub-Poisson because satellites are a minority of 
the galaxies. Lastly, when modeling the red and blue galaxies 
separately, P06 does not require that the number of central red 
galaxies and central blue galaxies be < 1 at fixed halo mass. 
P06 stipulate that if {N)m > 1, the 'first' galaxy is central and 
the remaining objects are satellites. However, this is done for 
both red and blue occupation functions. The implication is 
that there are two central galaxies per halo, one red and one 
blue, above Mo.red. When implementing an HOD that is sim- 
ilar to equations ([T]l and Q, P06 see no difference in the 
for the blue galaxies, but the f or the red galaxies increases 
by 20. See [Tinker et al.l (l2009al) for a full discussion of the 
implications of an HOD model of the type utilized by P06. 

4.5. Parameter Constraints for All the Samples 

Figure|2]shows the constraints in the Mred/Mbiue-/Rsat plane. 
The X-axis does not extend below Mred/A^Wue = 1 because that 
is the minimum value allowed from equation (O. For the UDS 
data, the lack of significant constraints is clear from the strong 
degeneracy between fusat and Mred/A^biue; with a larger frac- 
tion of red satellite galaxies, the need for a large ratio between 
the masses of red and blue central galaxies is reduced. This is 
a consequence of the small volume of the sample and the use 
of angular clustering. 

For the next three samples, the results are markedly differ- 
ent. All of the samples are consistent, within the l-cr con- 
straints, with red central galaxies being a random subsample 
of all central galaxies. A value of M^ed/Muue = 2 is excluded 
at more than 2-a for all the samples. The fraction of red satel- 
lite galaxies is high, being consistent with 60% in all sam- 
ples. Note that /Rsat = 0.6 is consistent with the results from 
the UDS data, in which case Mred/Mbiue ~ 1. Thus, if there is 
no strong evolution in the red satellite fraction from z ~ 0.9 to 
1 .4, the UDS data also suggest that red central galaxies are a 
random subsample of all galaxies. 

5. THE TIMESCALE FOR SATELLITE TRANSFORMATION 

5.1. A Simple Model 

Galaxies form within dark matter halos; the halo occupa- 
tion formalism is built upon this concept. High resolution 
collisionless N-body simulations track dark matter halos af- 
ter they have been accreted onto larger halos. Models that 
associate these subhalos with the satellite galaxies that con- 
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TABLE 1 
Parameters for Best-Fit Models 



Sample 


/Rsat 


/klTlLiX 


K 


logMsat 


log/3 


logMmin 


"red 


"blue 




V 




UDS 


0.22 


0.98 


2.47 


12.803 


-1.815 


11.838 


1.10 X 10" ^^ 


3.17 X 10" 


-3 


19 


10.8 


DEEP2 -20.5 


0.66 


0.25 


0.04 


13.365 


-0.363 


11.882 


8.46 X 10"^ 


2.13 X 10" 


-3 


29 


64.4 


DEEP2-19.5 


0.52 


0.25 


0.45 


12.800 


-1.891 


11.682 


2.37 X 10-3 


8.00 X 10" 


-3 


30 


63.7 


COMBO-17 


0.59 


0.17 


0.13 


12.694 


0.050 


11.391 


3.96 X 10-3 


1.23 X 10" 


-2 


46 


86.1 



Note. — The units of Mmin and Msat are h . v is the number of degrees of freedom in the sample. 
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Fig. 10. — Timescale for galaxies to transition from blue to red, or from 
star-forming to quiescent. All en'or bars are 2(j. The middle four (red) data 
points represent, in order of increasing redshift, COMBO-17, DEEP2 faint, 
DEEP2 bright, and the UDS. The blue data points are taken from other pa- 
pers. The datum at z = 2.3 is from the Itinker et al J <2009a) analysis of the 
clustering of DRGs. The upper datum at z = 0. 1 is from Wang et al. 1 2007), 
in which the exponential decay of star formation in satellites was constrained 
to be 2.5 Gyr We assume that satellites are "quenched" when the star forma- 
tion rate decreases by a factor of 10, or 2.3 e-folding times. The lower datum 
at z = 0. 1 is an estimate of the timescale for transitioning to the red sequence 
using the clustering data and analysis in Zehavi e t ah (2003). The green cir- 
cle at z = 1 .4 is the quenching timescale infen'ed from the UDS data if we 
assume that Mfej/Mbiue = in agreement with the other three samples. The 
shaded band is a power-law in which the quenching timescale varies with the 
dynamical time of dark matter halos, (1 -Fz)"'-^. 



stitute groups and clusters have been successful at reproduc- 
ing measured c lustering of galaxies at various redshifts (e.g., 
Kravtsovetanr2 004; Conrov et al. 2006; Wang et al. 200^ 
2007t iMarm et alj|2008t IMoster et al.ll2009t IWetzel & Whit3 



20091) . In these models, the known unknown is the exact re- 



lationship between halo and subhalo mass and the luminosity 
(or other galaxy properties, such as stellar mass) of the galaxy. 
A simple and successful approach has been to assume a mono- 
tonic relationship between galaxy luminosity and halo mass, 
where the exact functional form of the mass-to-light ratio is 
constrained to match the luminosity function of galaxies. Al- 
though there is kn own to b e some scatter in halo mass at fixed 
galaxy properties (iMore e t al. 2009), this monotonic approach 
has proved successful and useful in myriad studies, producing 
results quantitatively similar to those that incorporate scatter 
In the results in §4, we have made no assumptions about 



the population of dark matter subhalos within distinct halos. 
We have also made few assumptions about the MjL ratio as a 
function of halo mass because the samples we have analyzed 
have are defined with a threshold luminosity, rather than being 
broken into multiple luminosity bins. However, to make infer- 
ences about the evolution of satellite galaxies we must make 
use of N-body simulations that keep track of merger rates and 
the lifetimes of subhalos. For each observational sample, we 
have determined the total number of satellite galaxies in each 
halo, as well as the fraction of those galaxies that are clas- 
sified as red or quenched. To determine the quenching time 
of a galaxy once it is accreted, a simple approach is to con- 
sider all subhalos in a given halo that contain galaxies in the 
sample. Using the monotonic MjL ratio, this would simply 
be the (A^sat)M most massive subhalos (where the subhalos 
are ranked by their mass at the time of accretion, not their 
present-day mass, w hich is affected by orbit-dependent tidal 
stripping; s ee, e.g., IConrov et a l.1 l2006t i Moster et alj 120091; 
IWetzel & W hite 20oC These subhalos are then ranked by 
the time elapsed between accretion and the redshift of obser- 
vation. Assuming that the oldest subhalos are the ones that 
have had their star formation quenched and migrated to the 
red sequence, the value of /Rsat sets the quenching timescale. 

One complication is that some satellite galaxies may have 
been red before they became satellites. The results from all 
the samples in §3 indicate that, even at low halo masses, 
the fraction of field red galaxies is 20-30%. The mecha- 
nism by which these galaxies transitioned to the red sequence 
is unrelated to ram-pressure stripping or tidal stripping. A 
model in which field red galaxies are quenched by tidal in- 
teractions with nearby massive structures is also disfavored; 
an environmentally-dependent quenching mechanism of this 
type would enhance the large-scale clustering of red galax- 
ies. We have already noted that the relative bias of red to 
blue galaxies is difficult to match even if the probability of a 
central galaxy being red is independent of environment. An 
environmentally-dependent /Rcen is also disfavored at z = 
drinker et al. 2008b). 

To remove the contribution of satellite galaxies that were 
red before accretion, we determine the quenched fraction of 
satellite galaxies, /q, by 



/RsatHs; 



(8) 



*piev 



where «sat is the number density of satellites and Wprev is the 
number density of previously quenched satelhtes, given by 



= / dMn(M) 



f/Msub«sub(Msub|M)/Rcen(Msub) X (A^cen)Msub 

(9) 
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where nsub(A^sub|A/) is the subhalo mass function (where once 
again Msub is defined as the mass at the time of accretion), 
and /Rcen(A^sub) is the central red fraction determined for each 
sample. Note that nsub(A/sub I^C^ in our definition, does not 
have units of volume but rather is the number of subhalos 
within a given parent halo of mass M. We will describe our 
fitting function for this quantity in the following subsection. 
Here, we assume that the central galaxy red fraction does not 
strongly evolve since the time of satellite accretion, supported 
by the mild evolution in /Rcen across samples of §3. Equation 
^ yields a quenched fraction that is somewhat smaller than 
the overall red fraction. For three of the four samples, /rcbm 
is nearly independent of mass and Wprev ~ /Rmax^sat- For the 
best-fit model obtained for the UDS, /Rcen depends strongly 
on mass. However, in this model, the difference between f^^-dt 
and /q is minimal; in that particular model, /Rcen = below 
M ^ 10'^^ /!~'M0. Most satellites in this sample are below 
this mass, thus the fraction of subhalos that became red be- 
fore accretion is small. Recall, however, that the constraints 
on /Rcen in the UDS are broad. 

5.2. N-body Simulations 

To implement the model in the previous subsection, we re- 
quire two quantities: the subhalo mass function and the dis- 
tribution of accretion times for subhalos. We use the halo 
and subhalo catalo gs from the high - resolu tion N-body sim- 
ulation analyzed in IWetzel & Whit j (l2009h . The simulation 
is of the same cosmology as that assumed here. The vol- 
ume of the simulation is 200 /;"' Mpc per side, containing 
1500-^ particles, with a particle mass of 1.6 X IQ^Ii-^Mq and 
a force resolution of 3 /i"' kpc. Halos are identified using the 
friends-of-friends halo finding algorithm with a linking length 
of 0.168 times the mean interparticle separation. Halos de- 
fined in this way have dif ferent masses th at halos defined by 
A = 200 (White 2002; Ti nker et alJl2008ah . as done in all an- 
alytic calculations in this paper We con vert from FOF mass 
to A = 200 mass with a fac tor of 1.3. See lWetzel etal] (l2009h 
and lWetzel & Wliit3 (120091) for details of subhalo finding and 
tracking. 

Figure [8] shows the results for the subhalo mass function at 
six different redshifts between z = 0.9 and z= 1.6. The re- 
sults are plotted as a function of Msub/^host, in which units 
dN /dlnMsuh is independent of host halo mass. This self- 
similarity of the subhalo mass function has been shown be- 
fore for subhalo masses measured at a fixed redshift (e.g. , 
Gao et all 120041; iDe Lucia et all 12004 iKravtsov et all 120041; 
Wetzel et al]|2009l) . but note that Msub here is the mass at the 
time of accretion, the quantity that correlates with the ob- 
servational properties of the galaxy contained within them. 
Figure [8] also demonstrates that the subhalo mass function 
changes little with redshift. For each bin in Mhost, the N-body 
data turn over at low mass ratios. This is due to the finite 
particle number in the subhalo catalogs utilized; subhalos of 
the same initial mass will experience differential tidal forces. 
This causes some subhalos to drop below the particle limit of 
the subhalo catalog, thus the downturn is artificial. The region 
where the curves overlap is the true subhalo mass function. 
We use a fitting function of the form 



dN 



= 0.13 



sub 



-0.7 



c/lnMsub V^host 
when calculating equation (|9|l. 



exp 



-9.9 



/Ms. 



2.5 



(10) 
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Fig. 1 1. — The effect of imparting a critical mass scale of /Rcen for the 
faint red DEEP2 sample. The dotted curve is the model fit from Figure [5] 
The solid curve is the resulting coiTelation function from a model in which 
red galaxies rapidly dominate the central galaxy population above the critical 
mass obtained in the shock-heating model of Dekel & Birnboim (2006). The 
value of Mied/'Wbiiic is 10 in this model. The dashed curve is a less extreme 
model in which the majority of galaxies are red above the critical mass. For 
this model, -Mred/A^Wue = 5. 



Figure |9] shows the age distribution of subhalos at three dif- 
ferent redshifts. The y-axis is the cumulative distribution of 
satellites sorted by the time since they were accreted. The 
value on the y-axis indicates the fraction of subhalos that were 
accreted more than time t ago at each redshift. We define fsat 
as f(z)-f(accretion), and f(accretion) is defined as the time at 
which a subhalo is first linked to its parent halo by the FOF 
halo finder. The shape of P(fsat > t) clearly depends on red- 
shift; at earlier epochs, halo densities are higher by (1-l-z)-', 
thus the dynamical timescale decreases by (1-l-z)"'-^. Once 
the quenched fraction of satellites is determined, the quench- 
ing timescale,fQ found by smoothly interpolating between the 
numerical data in the P(fsat > distribution. We determine tq 
at the median redshift of each sample. In tests with the UDS 
sample, we find negligible differences between calculating tq 
at the median redshift and using the full N{z) of the sample. 

5.3. Results 

Figure [TOl shows the results of the model described in §5.2 
for all four samples (middle four points), with 2a error bars 
shown. All four samples are consistent with a value of tq ^ 
1.8 Gyr. There is no obvious dependence on redshift within 
these four samples. Figure [TO] also shows quenching times 
from highe r and lower redshift obtained from the literature. 
At z ~ 2. 3. [Tinker et al.l (|2009a) have analyzed the clustering 
of distant red galaxies (DRGs). These are massive galaxies 
that reside in halos with M^in ~ 10'^ Ii'^Mq. Some frac- 
tion of these galaxies have highly attenuated st ar formation 
rates, enough to be classified as red-and-dea d (iLabbe et alJ 
I2OOI iPanovich et aLlIlM iKriek et alJl2006h . iTinkeretalJ 
(I2009ah" found that the fraction of satellite galaxies that must 
be DRGs is between 0.5 and 1 {2a), while current estimates 
put t he fraction o f DRGs that are q uenched between 10-30% 
(iLabbe et al.ll2005; Papovich et al.ll2006.) . Using these errors 
and the model described above, the quenching time for these 
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galaxies is between 300 and 950 Myr (2a). 

At low redshift, current estimates for the quenching 
timescale of satellite galaxies are much longer IWang et all 
(|2007) used the clustering of star-forming and quiescent 
galaxies, as defined by the D„4000 break in the galaxy spec- 
trum (e.g., Kauffmann et al. 2003a ). to determine the quench- 
ing timescale of satellite galaxies in the SDSS. They found an 
e-folding time for the star formation rate of satellite galaxies 
of 2.5 Gyr Setting the quenching time to be the time at which 
a galaxy's star formation rate is attenuated by a factor of 10, 
we find fg = 5.75 Gyr This factor of 10 is roughly the differ- 
ence between the star-forming and quiescent samples in the 
UPS (W09), as well as the blue and red s ubsamples in DEEP2 
jNoeske et al.ll2007l) . Although lWang et al. (2006) provide no 
formal error estimate on this number, a rough error bar is 25 % . 

The blue circle at z = 0.1 is the quenchin g time we deter - 
mine through the clustering analysis of Zeh avi et al.l (l2005h . 
who use HOD modeling of color-selected samples in the 
SDSS to determine the red satellite fraction and the red cen- 
tral fraction of galaxies. For galaxies in magnitude bins cen- 
tered on M, =-19.5 and Mr = -20.5, the model in section §5.1 
yields quenching times of ^ 4 Gy r This est i mate is meant 
to be a consistency check on the iWang et al.l (l2006l) results; 
a more thorough analysis of the final data release of SDSS 
galaxies using these techniques will be presented by Wetzel 
et. al. (in preparation). The green circle at z = 1.4 is the 
quenching time obtained if we impose a prior of M^ed / A^biue < 
1.1 on the MCMC analysis of the UDS clustering. Such a 
prior increases the fraction of satellites that are red in order to 
match the clustering amplitude of red galaxies, thus decreas- 
ing the quenching time. 

The shaded band in Figure [TO] is a power-law that goes as 
(1 following the redshift dependence of the dynamical 

time of dark matter halos. From z = 0toz = 2.3, the results 
are consistent with a quenching timescale that is proportional 
both to the dark matter halo dynamical time and the satellite 
infaU time. 

6. DISCUSSION 

6.1. How do Central Galaxies Arrive on the Red Sequence ? 

The majority of galaxies in the universe are central galax- 
ies. As discussed in the introduction, current theories on the 
buildup of the red sequence focus on how to transform these 
central galaxies from star-forming objects to quiescent, red- 
and-dead objects. One set of ideas posits that there is a crit- 
ical halo mass at which star formation is quenched, either 
from AGN heating or through sustainable sh ock heating (eg, 
ICroton et al.1l2006HDekel & Birnboimll2006h . This theory of 
a critical halo mass is clearly at odds with the clustering re- 
sults analyzed here. If there exists a sharp critical mass scale 
above which galaxies become red, the relative clustering of 
red and blue galaxies would be significantly higher than that 
observed. 

Figure [TT] demonstrates this point quantitatively. The clus- 
tering data for the red galaxies in the faint DEEP2 sample are 
shown with the best-fit model from Figure|5] The solid curve 
is a model in which the fraction of red central galaxies transi- 
tions from to 1 sharply at M^n = 2 x 10'^ /;"' Mq, yielding 
A^ied/A^biue = 10. The shapc of /Rcen IS similar to the dashed 
curve from Figured To keep the number density fixed at the 
observed number density of red galaxies, f^^^x is reduced from 
0.6 to 0.45. But even with the lower satellite fraction in this 
new model, the amplitude of clustering is significantly higher 



than the data. The dashed curve shows the results for a less 
extreme model with a lower Mem and /Rmax = 0.75, yielding 
Med/A^biue = 5. This model is in better agreement with the 
data in the one-halo term, but is just as discrepant in the two- 
halo term. This type of model also fails for the bright DEEP2 
sample as well as the COMBO-17 data. For the UDS, a hard 
critical threshold is the best-fit model (cf. Figure [3]), but the 
errors are too large to distinguish between this model and one 
in which Mj-ed /Mbiue = 1 ■ 

The 'secular' scenario (Bo wer et al.ll2006l) . where instead 
of a critical halo mass threshold, there is a critical galaxy 
mass threshold, also predicts a strong dependence of /Rcen 
on halo mass. At M = 10"-^ /j-'Mq /Rcen 0, and by 
M = 10'^ ° Mq the f raction of red central galaxies is near 
unity (see Figure 6 in 'HopkinsetaLl l2008at) . The merger 
scenario has a weaker dependence of /Rcen on M, but still 
produces a high Mred/Mbiu e ratio. The halo merger rate de- 
pends httle on halo mass (ICohn et al.ll2"00H iFakhouri & Mai 
l2007h . This is also true of the merger rate of subhalos with 
host halos, which should more accurately track the galaxy 
merger rate ( Wetzel et al. 2009). However, the mapping be- 
tween halo mass and galaxy mass is complex; the Mgai/M ra- 
tio peaks at M ~ 10^^ /z~'M0 and falls off as a power law 
at higher and lower halo halo masses (see, eg. |Wang et alj 
l2007tlConroy & Wechsled2009HMoster et al.f 2009). Because 
Mgai ^ M°-4 at M > 10'2 h-^M.Q, a 10:1 halo-halo minor 
merger can contain gala xies with mass ratios of 3:1, yieldin g 
a major galaxy merger (iHopkins et al.ll2008aL iMalleilllOOSh . 
In all calculations bel ow we use the redshift-d ependent stel- 
lar mass functions of Marche sini et al] ("2007) assign stellar 
masses to halos and subhalos using abundance matching, in- 
terpolating between redshifts to obtain the stellar mass func- 
tion at the redshift of a given sample. 

Figure [T2l shows the z = 1.5 red central fraction induced 
by major mergers, which we define as stellar mass ratios of 
3:1 or closer, that have occured over the previous Gyr The 
small circles indicate the fraction of halos that have expe- 
rienced a halo-subhalo (central-satellite) merger with mass 
ratios smaller t han 2:1, 3:1 , and 5 :1, respectively, from the 
simulations of IWetzel et alJ (l2009h . The lack of strong de- 
pendence on halo mass is clear Using abundance matching 
to map halo mass onto galaxy mass, a 3:1 mass ratio for 
galaxies in those same halos yields a red central fraction of 
- 10% at M = 10'2 /z"' Mq, rising to - 30% at lO'^ /i"' Mq 
and - 60% at 10"* /z"' M©. The right hand panel of Figure 
fT2l compares the central occupation function for red galax- 
ies for model with Mred/Mbiue = 1 to (A^cen)M created by the 
merger scenario. At z = 1 .5, the merger scenario yields a cen- 
tral galaxy number density that is too low by roughly a factor 
of 3. As stated above, the halo merger rates are calculated 
using a lookback time of 1 Gyr Increasing this timescale in- 
creases the fraction of halos that have experienced mergers, 
but it does not change the shape of the merged fraction with 
halo mass, thus we consider the amplitude of this function to 
be a free parameter set by the overall number of red galax- 
ies. Using the 1 Gyr lookback time to calculate /merge, the 
merger scenario yields M,-ed/Mbiue = 2.2. Increasing the nor- 
malization of /merge to producc the correct number density of 
yields Mred/Mbiue = 2.9. These values are in agreement with 
the weak constraints on /Rcen shown in Figure [T] At z = 0.8, 
the redshift of the faint DEEP2 sample, the merger scenario 
yields Mred/Mbiue = 3.8 using the 1 Gyr lookback time, and it 
yields Mi-ed/Mbiue = 5.2 if the amplitude of /merge is increased 
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Fig. 12. — Left Panel: Small circles with errors represent the fraction of halos that have experienced a halo-subhalo (central-satellite) merger within the last 
Gyr with mass ratios smaller than 5:1 (upper), 3:1 (middle), and 2:1 (lower). The dotted lines show the mean values as for each mass ratio limit. The thick solid 
curve shows the fraction of halos whose ce ntral galaxy h as experienced a merger with a stellar mass ratio smaller than 3:1. The large red circles represent the 
fraction of red central galaxies from Zehav i^t al.l j2005l) analysis of SDSS galaxy clustering. Right Panel: Central occupation functions. The thin solid curve 
represents (Aten)M for the full UDS sample. The thick red line represents a model with Mred/^blue = 1- The thick green line represents (Ween)*/ produced by 
the merger scenario. 



to produce the proper red central number density. The higher 
values reflect the change in the halo mass function; under 
the merger scenario most high-mass halos are quenched, and 
the abundance of high-mass halos has increased substantively 
from z = 1.4 to z = 0.8. Results for COMBO-17 are similar. 
These values are strongly excluded for both the DEEP2 and 
COMBO-17 results. Thus, while the UDS data are broadly 
consistent with the merger scenario, DEEP2 and COMBO-17 
are not. 

To connect with the low redshift universe. Figure [12] also 
shows the depende nce of ./keen o n halo mass obtained from 
SDSS clustering by IZehavi et al.l (|2005). The shape of /R^en 
at z = 0.1 i s consistent with the merger scenario (see also 
iHopkins et a l. 2008a). If the COMBO-17 results are repre- 
sentative of the clustering of blue and red galaxies at z = 0.6, 
what occurs from z = 0.6 to z = to induce the change in /Rcen? 
Selection effects seem unlikely to impart a mass dependence 
in /Rcen in SDSS; indeed, the results of M ailer etal.1 (l2009h 
demonstrate that many low-luminosity galaxies in SDSS are 
misclassified as red due to dust effects. Similarly, selection 
biases seem unlikely to remove a mass dependence in /Rcen in 
COMBO-17 and DEEP2. Red galaxies are selected by differ- 
ent methods in these two samples; for COMBO-17, spectral 
energy distribution fitting is used; in DEEP2, a color cut sim- 
ilar to SDSS is employed. 

Is there evolution in the mechanism by which central galax- 
ies quench their star formation? For halos of < 10^'^ M0, 
the red central fraction for the high-redshift samples is simi- 
lar to z = estimates, modulo the differences in the classifica- 
tion of red galaxies between the various surveys. But /Rcen at 
M ~ io'^~ /;"' Mq must increase by nearly a factor of two, 
and this fraction must increase by nearly a factor of three at 
M> hr^MQ. More time elapses from z = 0.6 to z= 0(5.8 



Gyr) than elapses from z = 1 .4 to z = 0.6 (3.3 Gyr). Although 
halo dynamical times evolve strongly with redshift, there are 
roughly the same number of dynamical times between SDSS 
and COMBO-17 as there are between COMBO-17 and the 
UDS (approximately 6). Thus evolution should not be ruled 
out as a possibility. 

6.2. The Role of Satellite Galaxies in Building up the Red 
Sequence 

At all redshifts, from z = 2.3 to z = 0, the fraction of satel- 
lites that are red is roughly consistent with ^ 60% of all satel- 
lites in the sample. This is more than double the fraction of 
centrals which are red, demonstrating that satellites are not 
simply made red before infall. Figure [13] show what frac- 
tion of the red sequence is made up satellite galaxies for each 
sample, as well as the fraction of those satellites that were 
quenched after becoming satellites. At redshifts less than 
z = 2, 60% of satellites being red translates to ^ 30% of the 
red sequence being made of up satellites. Roughly 70-80% of 
red satellites became red after accretion, in agreement with 
results for low-mass galaxies at z = (Ivan den Bosch et aP 
i2008t) . although the UDS results are consistent with a value 
as low of 10%. Thus satellite galaxies play a major role in 
building up the red sequence. The surprising fact is that these 
numbers do not vary substantially with redshift. The survival 
time of a subhalo scales with the dynamical time of the host 
halo, (l+z)"'^. Thus, to keep /Rsat a constant requires the 
same redshift dependence in the quenching time, as shown 
in Figure [TO] Given the finite lifetimes of subhalos, satel- 
lite galaxies observed at various redshifts are distinct sets of 
galaxies. The fraction of subhalos at z = 0.6 that survive to 
z = is ~ 20%. The fraction from z = 0.9 that are extant to- 
day is less than 5%. Thus the satellite galaxies analyzed in 
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the SDSS by I van den Bosch et all (l2008h are not the same as 
th ose quantified here , yet their statistics are quite similar 

IWang et aTl (l2007h used clustering to constrain the quench- 
ing time at z = (cf, FigurefTOli. They posit that this quenching 
time does not evolve with redshift, and that this fixed tq value 
yields a color-density relation that largely goes away at z > 1, 
in qua litative agreement with the observations of lCooper et alj 
(l2007h . However, a quenching time that is fixed at 5 Gyr 
would be incompatible with all of the correlation functions 
analyzed in this paper. From Figure |9] it can be seen that no 
satellites are older than 6 Gyr at z > 0.6, thus there would 
be no satellites quenched after infall at z > 0.6. The satel- 
lite red fraction would simply trace the central red fraction, 
and the cluste ring of red g alaxies would be the same as blue 
galaxies. The Coope r et al.l (12007ft measurements were made 
using the DEEP2 survey, t hus they rn ust be consistent with 
the correlation functions of ICoil et al.l 12008) from the same 
data. The lack of a strong color-density relation, as well as 
its gradual attenuation with increasing redshift, are compati- 
ble with the HOD results obtained here; if red central galaxies 
are a random sample of all central galaxies, then they produce 
no color-density relation. The high fraction of satellite galax- 
ies that are red would induce a correlation between color and 
density, but as redshift increases the overall fraction of galax- 
ies that are satellites decreases. Thus /Rsat can remain fixed in 
redshift and still result in a decreasing satellite fraction within 
the red sequence. It should be noted, however, that the high 
relative bias between red and blue galaxies at 1 < z < 2 in 
W09 impli es that the colo r-density relation does exist in the 
UDS. Oua dri et al.1 (l2007l) find the same result at higher red- 
shift in the MUSYC survey. Further data are required to re- 
solve this discrepancy. 

The simple model for obtaining the satellite quenching time 
in §5.1 neglects certain effects. The quenching time may de- 
pend on both the age of a satellite and the details of its orbit (as 
suggested by Balogh et al. 2009), as well as the mass of the 
subhalo or the stellar mass of the satellite galaxy. Another ex- 
plicit assumption in this model is that the quenching can only 
begin after a galaxy passes inside the virial radius of the par- 
ent halo. These questions will be addressed in a subsequent 
paper (Wetzel et. al., in preparation). 

6.3. Inferences from the Luminosity Function 

At redshift zero, our now-standard picture is that the frac- 
tion of central galaxies that are red increases with halo mass. 
Above M ft! 10'^^ Mfr^, nearly all central galaxies are re d 
dZehavi et al.ll2005l:IWeinmann et alj|200a lYang et alj|2008l) . 
This scenario is influenced by measurements of the galaxy 
luminosity function; red fraction and galaxy l uminosity are 
strongly correlated (e.g.. iBlanton et aDl2003 lj others). The 
results at higher redshift, however, are quite different. Above 
Mb = -19.5, the abundance o f red galaxies in DEEP2 has little 
correlation with luminosity (IWillmer et al.ll2006) . In Figure 
[T] the red fraction increases by only 9% between the faint and 
bright DEEP2 samples. At lower luminosities, the fraction 
of red galaxies does decrease. Thus, the value of Mred/Mbiue 
should increase at lower luminosity thresholds. If there is a 
preferential halo mass scale for star formation quenching in 
central galaxies, it is below the mass scale probed by the sam- 
ples we have analyzed, M < 10" M0. 

6.4. Possible Systematic Errors 

A concern in this analysis is the use of diagonal error bars 
only, as opposed to using the full covariance matrix for each 
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Fig. 13. — Circles indicate the fraction of red satellite galaxies that were 
quenched after accretion onto their parent halo, /q, as defined by equation 
(Sj. Squares indicate the fraction of the red sequence made up of satellite 
galaxies at the magnitude threshold for each sample. The redshifts for the 
squares ai'e slightly offset for clarity. All en'orbars are 2a. 



sample. The clustering statistics used here are known to 
be correlated in the two-halo term and l argely uncorrected 
in the one-halo term (see, eg, figures in IZehavi et al J 120051 ; 
iPhleps etal.1 llooa lBlakeetal.1 l2008t iTinker et all I2009al) . 
Thus, the use of diagonal errors over-weights the contribution 
of the data in the two-halo term to the of each fit. Us- 
ing our 720 Mpc simulation, we have made mock galaxy 
distributions that roughly match the faint DEEP2 sample. The 
mass resolution in this simulation is not enough to model M < 
10"-^ r'Mo halos, thus we are forced to modify (A^cen)M 
of the best-fit model to obtain the proper number density of 
galaxies. We divide the simulation volume into 6-^ = 216 cu- 
bic subvolumes of 120 Mpc per side. We calculate Wp(ry,) 
within each subvolume using the Landy-Szalay estimator and 
accounting for the variation in the number density of galaxies 
in each subvolume. We are not able to obtain a robust co- 
variance matrix for both the red and blue subsamples, thus we 
adopt the same (normalized) covariance matrix for both sub- 
samples. We multiply each element by the diagonal elements 
measured by ICoil et af] (l2008h . ie C,y = C/jg/e;, where Cij is 

the normalized matrix from the mocks, Cij = Cijj ^CuCjj. Al- 
though the covariance matrices of the red and blue subsam- 
ples should vary in detail, our test replicates the general con- 
sequences of a covariance matrix: some correlation of data 
in the one-halo term, with an increasing amplitude of off- 
diagonal terms in the two-halo term. We find virtually no dif- 
ference in the Mred/Mbiue-/Rsat Constraints when using these 
test matrices. We conclude that our exact answers may vary 
slightly with use of proper estimates of Cy , but the results will 
not change enough to alter our conclusions. 

Selection effects are also a concern. If /Rcen did depend on 
halo mass, either through a critical mass scale or a weaker 
dependence from a merger-like scenario, dust contamination 
would tend to mask this mass dependence. If galaxies are mis- 
classified as red or quenched due to dust contamination, the 
large-scale amplitude of red galaxy clustering is attenuated. 
In terms of the HOD models presented here, and under the 
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assumption that /Rcen was actually monotonically increasing 
with halo mass, dust would take low-mass blue central galax- 
ies and place then in the red sample, reducing Mred/A^Wue- 
However, this effect seems unlikely to produce /Rcen in the 
best-fit DEEP2 model. First, if /Rcen had a strong mass de- 
pendence to it, such that the fraction of red central galaxies at 
M > 10'-' /i"' Mq were > 75% (a value taken from the anal - 
ysis of z = clustering in the SDSS by IZehavi etal1l2005h . 
then a mass-independent dust contamination would smooth 
our /Rcen but not be able to remove this trend completely. 
More importantly, the low values of /Rcen at high halo masses 
in the DEEP2 models 0.15-0.25) are at odds with this 
scenario. 

One degree of freedom we have not introduced into the 
HOD model used here is a mass-dependence on the fraction of 
satellites that are red. We have performed tests with the faint 
DEEP2 data, introducing a mass scale below which /Rsat is at- 
tenuated with a Gaussian cutoff. The constraints on this cut- 
off mass are weak and do not change either tp or M red/A^biue- 
The color-dependent HOD in IZehavi et alJ (l2005i) contained 
a mass-dependent parameterization in /Rsat, but in most clus- 
tering samples /Rsat was found to be constant. If red satellite 
galaxies prefer higher mass halos, the red correlation func- 
tion would contain a more pronounced break between the 1 - 
halo and 2-halo terms, and the shape of the red correlation 
function at small scales would not be a power-law. In the 
DEEP2 data, Wp(rp) for the red galaxies is a power law down 
to rp ^ 50 kpc, and it does not contain a strong break at 
rp^ I /;"' Mpc. 

Although Qfsat = 1 is a well-motivated prior for luminosity 
thresholds, assuming this power-law index for color subsam- 
ples is less cut and dry. A better fit to the faint DEEP2 data 
is obtained if I parameterize (A^sat)M for blue and red subsam- 
ples as independent power-law functions, with different val- 
ues of Msat and a^^t for each. The for this model is 38.6, 
a significant improvement on the fit from the fiducial model. 
However, to produce this x^, both asat for the red and blue 
galaxies is well below unity. If asat = 1 for the full sample, 
then a value of Ofsat < 1 for blue galaxies requires cksat > 1 for 
the red subsample (or vice versa). Models that produce a bet- 
ter fit than the fiducial model have both values below unity, 
with the best-fit model yielding Ofsat = 0.39 for the red sub- 
sample and Q!sat = 0.48 for the blue subsample. These values 
are firmly excluded from theoretical expectations and obser- 
vational data (cf. §3.1). Within the range 0.9 < Ofsat < 11, 
our constraints in Figure|7]are unchanged. 

6.5. Comparison to Results from the NDWFS 

iBrown et al.l (l2008l) perform HOD analysis on the cluster- 
ing of red g alaxies from z = 0.2- 1.0 in the NDWFS imag- 
mg survey (lJannuzi&Deylll999t) . At - 7 square degrees, 
this sample is much lar ger than all th e samples analyzed in 
this paper Additionally. iBrown et alj (|2008) use the full co- 
variance for all halo occupation analysis. They find that the 
/keen = 0.5 at lO"'^ M0, increasing with halo mass. How- 
ever, the HOD model for central galaxies used in their anal- 
ysis is identical to Equation ([T]) here. Essentially, they model 
red galaxy samples with an HOD constructed for a full galaxy 
sample. This means that f^cen must increase with mass and 
it will approach unity at high halos masses by construction. 
Without a corresponding sample of blue galaxies, it is difficult 
to const rain the asyrnntotic fraction of red central galaxies. 
Because [Brown etaTI (l2008l) only have photometric redshifts 



of the galaxies in their sample, they can only measure w(9), 
which has larger fractional errors at large scales relative to 
measurements of Wp{rp) over similar volumes. All HOD fits 
to their data are somewhat high in the two-halo term, which 
could be caused by the assumption that /Rcen increases to unity 
(cf. Figure fTTI). However, the errors on w{6) are large enough 
such that they are all statistically good fits. Corresponding 
measurements of the clustering of blue galaxies in the same 
fi eld would help resol ve the discrepancy between the results 
of lBrown et alj ( l2008l) and those in this paper. 

7. SUMMARY 

In this paper we have analyzed clustering measurements 
of red and blue galaxies at four different redshifts, z= 1 .4, 
z = 0.9, z = 0.8, and z = 0.6. We have used the halo occu- 
pation distribution to determine the fraction galaxies are that 
red, broken into contributions of red central galaxies and red 
satellite galaxies. In all the samples analyzed, the higher clus- 
tering amplitude of red galaxies can be fully accounted for 
the the high fraction of satellite galaxies that are red; on aver- 
age, the overall red fraction of each sample is ^ 25%, while 
the fraction of satellites that are red is ~ 60%. Thus the frac- 
tion of the red sequence that is satellite galaxies is ~ 30%. 
The high value of /Rsat produces enhanced clustering at both 
large and small scales relative to blue galaxies above the same 
magnitude threshold. Imparting a halo mass dependence on 
the fraction of central galaxies that are red enhances the rel- 
ative bias between red and blue subsamples are large scales. 
The data are not compatible with such a scenario, favoring 
Mred/M\,iue = 1 for Central galaxies. 

From galaxy-halo mass-to-light ratios (and stellar mass 
to halo mass ratios), it is inferred that galaxy formation 
efficie n cy i s maximal at M ~ 10 '^ /i"' M^ (Yang et alj 
200l Ivan den Bosch et aP 120071: [T inker etal. 2001 
Conrov & Wechslen 12009^ iMoster et al.l i2009i) . This 
mass scale increases slowly with red shift out to z ~ 1 
(IConroy & Wechsled[2009[: iMoster etan '2009). This implies 
that a central galaxy's integrated star formation history is 
dependent on halo mass. This leads to a correlation between 
instantaneous star formation rate and halo mass (see Figure 8 
in Conroy & Wechsler 2009). In contrast, our results imply 
that the fraction of galaxies that are quenched is independent 
of halo mass at z > 0.6. 

Models in which central galaxies become red at a critical 
halo mass scale of ^ 10'^ h~^MQ are strongly excluded by 
the data available at z < 1 . A scenario in which mergers pro- 
duce red galaxies has a weaker dependence of red central frac- 
tion on halo mass, but these models also produce values of 
Med/^biue larger than unity, values that outside the 2-a con- 
straints from DEEP2 and COMBO- 17. 

The constant value of /Rsat implies that the quenching 
timescale for satellite galaxies depends on redshift in the 
same manner as the dynamical timescale of dark matter ha- 
los. Given the errors, tq for all four samples is consistent with 
fQ = 1 .8 Gyr. Current estimates of fg at z = are much longer, 
~ 5 Gyr, while preliminary estimates of Jq at z > 2 are less 
than a Gyr. 

Given the small volume of each sample, more data are re- 
quired to ameliorate the possibility that sample variance is 
significantly influencing our results. Additionally, clustering 
measurements at redshift between and 0.6 can shed light on 
discrepancies seen between z = and the samples analyzed 
here. Upcoming results from the G alaxy and Mass Assembly 
survey (GAMA; [Driver et"al]|2009l) . as well as the clustering 
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of blue galaxies from the NDWFS, should show intermedi- 
ate stages of evolution in both the quenching mechanism for 
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